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摘要 利用 星系 解构 软件 GALFIT 通 过 面 亮 度 轮廓 拟 合 对 近邻 早 型 旋涡 星系 M81 (NGC 3031) 进 行 形态 学 解 
TT FJ, 旨 在 探究 M81 星 系 的 结构 组 成 并 对 其 进行 形态 学 量化 . 通过 6 种 解构 模式 , 对 M81 进 行 了 不 同 复杂 程度 的 结 
构 分 解 , 其 中 最 复杂 的 解构 模式 包含 核 球 、 盘 、 外 旋 臂 、 内 旋 臂 、 星 系 核 5 个 子 结构 . 研究 结果 显示 , M81 
一 个 Sérsic 指 数 约 为 5.0 的 经 典 核 球 , 其 形态 和 光度 在 不 同 解构 模式 中 均 保 持 稳 定 ; M81 星 系 盘 的 S6rsic 指 数 约 
231.2, 但 它 的 形态 参数 和 光度 与 是 否 分 解 内 旋 辟 相关 . 不 同 子 结构 的 组 合 对 作为 混合 体 的 星系 整体 的 形态 有 不 
可 忽视 的 影响 . 星系 解构 的 结果 提供 了 不 同 解构 模式 适用 性 的 建议 : 其 中 核 球 + 盘 + 星系 核 的 三 成 分 解构 适用 
于 大 样本 星系 的 核 - 盘 研 究 ; 而 考虑 旋 臂 的 复杂 解构 则 适合 于 对 星系 子 结构 的 精确 测量 , 如 小 样本 (或 个 源 ) 研 究 . 
基于 Spitzer-The Infrared Array Camera (IRAC) 4.5 hm 的 单 波 段 图 像 的 形态 学 解构 研究 是 后 续 一 系列 研究 的 
开始 , 在 此 基础 上 未 来 将 会 对 M81 进 行 多 波段 解构 , 同时 研究 不 同 子 结构 的 光谱 能 量 分 布 和 星 族 性 质 , 并 推断 
M81 各 子 结构 的 形成 历史 和 演化 过 程 . 
关键 词 星系 : 个 别 : M81, 核 球 , 星系 : 基本 参数 : 形态 学 , 星系 : 结构 , 星系 : HR 
中 图 分 类 号 : P157; ”文献 标识 码 : A 


1 引言 过 程 T 在 伴星 系 或 
星系 的 形态 是 星系 最 直接 的 观测 表征 . 决定 ”其 他 不对 称 引力 势 的 影响 F, 可 能 产生 旋 辟 或 其 
星系 形态 的 主要 因素 是 星系 的 形成 和 演化 ("31. 例 他 子 结 构 , 恒星 运 Bo ia ci 速度 弥散 
C 如 椭圆 星系 (或 球状 星系 ) 通 常 被 认为 是 通过 短 时 DEO. 不 同 的 形成 过 程 在 星系 中 产生 了 不 同性 质 
标的 剧烈 过 程 (如 并 合 ) 形 成 的 , 这 种 合并 的 机 制 加 KEE. 椭圆 星系 中 主要 是 年 老 恒星 族 , 恒星 形成 
XE T BL £ BA RS BRR Pet ”活动 水 平 低 , 星系 整体 颜色 偏 红 ; 盘 星系 中 的 恒星 
E, 并 且 打 破 原本 星系 内 部 秩序 的 恒星 轨道 从 而 形 。”” 族 通常 更 年 轻 , 星系 整体 颜色 偏 蓝 , 恒星 形成 率 相 
成 球体 结构 , 整体 形态 相对 光滑 , 没有 断裂 或 其 他 ”对 更 高 . 研究 星系 的 形态 是 认识 星系 中 恒星 族 性 
子 结构 ; 恒星 运动 以 随机 运动 为 主 , 速度 弥散 比较 i, 了 解 星系 形成 和 演化 的 基础 . 
AUD, 盘 主 导 的 星系 则 被 认为 是 由 长 时 标的 吸 积 一 个 星系 也 常常 由 多 个 不 同形 态 的 子 结构 组 
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成 , 尤其 是 对 于 盘 星系 来 说 , 其 主要 子 结构 是 核 球 、 其 中 ,表示 系统 中 心 的 表面 亮度 , rs 表示 盘 的 标 
盘 和 旋 臂 . 盘 星 系 的 核 球 有 着 和 椭圆 星系 相似 的 度 长 度 (Scale Length). 

观测 特征 , 它们 有 相似 的 面 亮度 轮廓 Po 9; 这 意味 然而 , 无 论 是 用 r1A4 定 律 拟 合 球状 系统 , 还 是 用 
着 核 球 可 能 和 椭圆 星系 一 样 是 由 短 时 标的 剧烈 过 。 ”指数 定律 拟 合 盘 状 系统 , 通常 都 很 难得 到 理想 的 拟 
程 (如 并 合 ) 产 生 . 早 型 旋涡 星系 在 其 演化 早期 就 形 RAR. 在 这 种 情况 下 , 人 们 发 展 了 一 个 更 加 灵活 
成 了 核 球 , 而 盘 结构 可 能 是 通过 后 期 吸 积 气 体 而 产 的 函数 一 Sérsic 函 数 一 来 拟 合 星 系 的 面 亮度 的 , K 


生 的 子 结构 09. 旋 臂 结构 通常 是 星系 中 驱动 气体 的 数 形式 如 下 
重要 动力 来 源 , 旋 辟 会 聚集 星系 中 的 气体 , 从 而 会 pO <5 ov [^^ (^ ü J) (3) 
产生 大 量 的 恒星 形成 区 和 年 轻 恒星 3. 不 过 旋 臂 的 Te | 
成 有 不 少 争议 , 不 少 天 文学 家 认为 旋 臂 是 由 于 盘 其 中 指数 mn 称 为 <96rsic 指 数 "， 表示 亮度 的 集中 程 
受到 局 部 的 扰动 而 产生 的 39， RE; k 是 与 n 有 关 的 系数 (因此 不 属于 独立 的 自由 参 
对 于 由 不 同 子 结构 组 合 而 成 的 星系 , 其 整体 性 #0. 可 以 看 出 , Sérsic 函 数 把 r1/4 定 律 和 指数 定律 
质 是 不 同 子 结构 性 质 的 综合 , 也 就 是 说 , 星系 中 的 。 ”做 了 一 个 统一 也 就 是 当 Sersic 指 数 n 等 于 4 (此 时 


子 结构 会 影响 作为 各 种 子 结构 混合 体 的 星系 整体 的 < 等 于 7. 67) 时 ， Skrsic 函 数 就 是 rl4 定 律 ; Sérsic 
的 性 质 . 例如 , 一 部 分 绿 谷 星系 其 实 就 是 由 恒星 形 。 ”指数 n 等 于 1 时 ，S6rsic 函 数 就 是 指数 定律 (r。= 
成 率 高 的 盘 和 恒星 形成 率 低 的 核 球 平均 作用 的 结 1.678rs), Sérsic 函 数 把 指数 n 作 为 自由 参数 , 能 

RU 而 常见 的 星系 颜色 的 双 峰 性 也 是 由 于 核 球 和 更 加 灵活 地 拟 合 星系 面 亮度 的 径 向 轮廓 . 由 于 这 
盘 的 各 自 特征 的 表现 ””. 不 同 子 结构 的 混合 而 造成 — 样 的 灵活 性 , 因此 能 够 得 到 更 加 理想 的 拟 合 结果 . 
的 星系 物理 参数 估计 偏差 可 能 导致 对 星系 形成 和 目前 , S6rsic 函 数 是 最 常用 的 星系 面 亮 度 径 向 轮廓 
演化 理解 的 错误 . 如 果 对 星系 进行 解构 , 分 别 研究 。 拟 合 函数 , 包括 对 于 星系 不 同 子 结构 的 拟 合 . 星 
不 同 子 结构 各 自 的 性 质 则 可 以 有 效 地 降低 子 结构 。 系 形态 学 解构 软件 的 问世 使 得 人 们 能 够 非常 便利 


混合 带 来 的 不 确定 性 . 地 分 离 星 系 中 的 各 结构 成 分 并 对 其 进行 形态 学 量 
二 星系 子 结构 的 形态 进行 量化 是 星系 解构 的 。 ”化 [25-29]. 
基础 , 而 面 亮 度 径 向 轮廓 拟 合 则 是 最 常用 的 量化 子 然而 , 星系 形态 学 解构 结果 中 的 不 确定 性 是 该 


结构 形态 的 工具 , 因为 不 同形 态 的 恒星 系统 , 其 面 。 ”类 工作 的 主要 问题 之 一 , 解构 过 程 的 主观 性 和 可 变 
亮度 径 向 轮廓 遵 循 不 同 的 函数 形式 W**. 对 于 内 性 比较 强 . 不 同 目标 的 星系 形态 学 解构 其 模式 和 结 
部 恒星 以 随机 运动 为 主 的 球状 系统 (如 椭圆 星系 或 果 可 以 是 完全 不 同 的 . 对 于 盘 / 旋 涡 星系 , 人 们 通常 
HE ARIA HER, 它们 遵循 基本 得 理论 所), 其 面 ”采用 核 球 + 盘 的 解构 模式 , 不 过 这 种 两 成 分 模式 可 
亮度 径 向 轮廓 大 致 符合 de VaucouleurssE f£ (Bl r!/^ 能 过 于 简单 , 在 分 离 后 的 剩余 图 像 中 会 留 有 旋 臂 
~ ner) Pl, 函数 形式 如 下 : 等 其 他 子 结构 成 分 . ; 这 些 剩余 的 子 结构 如 果 没 有 
pai 被 纳入 独立 的 形态 学 成 分 是 否 会 对 拟 合 造成 比较 

us zap] id (=) i] (O0) 严 重 的 影响 ” 此外, 有 些 星系 的 中 心 有 活 动 星系 核 

其 中 , "为 星系 中 的 径 向 距离, 7 为 有 效 半径 , U(r) Gotive Galactic Nucleus, ACT 人 CN 人 和 会 

表示 在 径 向 距离 r 处 的 面 亮 度 ， VAR TEA BEE 何 影 影响 星系 的 形态 学 解构 结果 ; 同时 ， 不 同 的 解构 

v. Ab RYT SEO 模式 又 会 对 作为 独立 成 分 的 AGN 产 生 怎 样 的 影响 ? 

得 状 系统 的 面 亮度 径 向 轮 序 则 更 倾向 于 不 同 科研 目标 的 工作 应 该 要 如 何 选择 适合 自身 的 


遵从 指数 函数 , 因此 将 这 种 盘 结 构 称 为 “指数 ERRA? 


H 


fi (Exponential Disk)" 23, 函数 形式 如 下 : 为 了 回答 这 些 问题 , 我 们 开展 了 这 项 对 早 型 旋 
涡 星系 M81 (NGC 3031) 的 形态 学 解构 工作 , 以 量 
i c d Ge ij (2) 化 M81 的 形态 学 结构 组 成 , 并 探究 不 同 的 形态 学 解 
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构 模 式 对 星系 子 结构 拟 合 结果 的 影响 . M81 的 形态 
学 分 类 是 SAab 型 , 距离 地 球 约 3.6 Mpc (JKP ŽE), 
视角 尺度 约 为 26' x 14'. 近 距 离 、 大 角 尺 度 使 得 对 
M81 的 研究 可 以 在 本 征 高 空间 分 辨 率 的 条 件 下 展 
JF. 迄今 为 止 , 已 经 有 许多 工作 对 M81 的 多 种 性 质 
进行 了 研究 , 包括 星 族 成 分 B0 3、 人 尘埃 分 布 B3-361、 
S 态 学 B739 等 . 本 文 星 现 的 工作 是 对 M81 单 波段 形 
态 学 进行 的 深入 研究 . 以 这 项 工作 为 基础 , 我 们 后 
续 的 工作 将 对 M81 进 行 多 波段 形态 学 解构 , 并 系统 
研究 M81 各 子 结构 的 光谱 能 量 分 布 (Spectral Ener- 
gy Distribution, SED)、 星 族 性 质 以 及 星系 整体 和 
各 子 结构 的 形成 和 演化 过 程 . 

我 们 将 在 第 2 节 介绍 数据 来 源 和 图 像 处 理 的 方 
法 , 第 3 节 中 介绍 星系 形态 学 解构 的 过 程 , 解构 的 结 
果 将 会 在 第 4 节 展示 , 并 在 第 5 节 进 行 讨论 , 第 6 节 将 
对 本 工作 的 结论 进行 总 结 . 在 本 文中 , 1 的 物理 尺 
度 约 等 于 1 kpc. 


2 数据 来 源 及 其 初步 处 理 
2.1 ”数据 的 获取 

由 于 M81 是 著名 的 近邻 星系 , 因此 各 波段 的 观 
测 数据 都 很 丰富 . 在 本 工作 中 , 我 们 选择 分 辨 率 高 
受到 恒星 形成 区 、 不 规则 团 块 结构 、 人 尘埃 辐射 
影响 比较 小 的 Spitzer-The Infrared Array Camera 
(IRAC) 4.5 hm 的 图 像 数据 !. 数据 的 来 源 是 The 
Spitzer Survey of Stellar Structure in Galaxies 
(S4G) 项 H 4942, nf fENational Aeronautics and 
Space Administration (NASA)/Infrared Process- 
ing and Analysis Center (IPAC) Infrared Science 
Archive (IRSA) È P FE. 

KR f Spitzer-IRAC 4.5 pm 的 图 像 数 据 , 本 工作 
也 用 到 来 自 Galaxy Evolution Explorer (GALEX) 
Far ultraviolet (FUV) 的 图 像 数 据 , 作为 构建 旋 辟 
模型 的 参考 后]. 因为 旋 臂 主要 由 年 轻 恒星 族 组 成 ， 
选择 GALEX FUV 波 段 能 够 最 大 化 消除 年 老 恒 星 
族 对 旋 臂 结构 建 模 的 影响 . 该 观测 波段 数据 来 源 于 


che 


1Spitzer-IRAC 3.6 hm 的 M81 图 像 受到 饱和 像素 的 影响 , 因此 本 了 


*https://irsa.ipac.caltech.edu/frontpage/. 
Shttp://www.galex.caltech.edu/researcher/data.html. 
“https: //asd.gsfc.nasa.gov/archive/galex/. 
5http://github.com/kbarbary/extinction. 


GALEX data release 6 (DR6)?, nf 7EGALEX 
Guest Investigator? B |j rp F aX. 

在 本 工作 中 , GALEX FUV 图 像 只 用 来 作为 旋 
臂 参数 建 模 的 参考 , 星系 的 全 部 形态 学 解构 过 程 以 
及 所 有 结果 都 是 基于 Spitzer-IRAC 4.5 hm 图 像 . 图 
1 是 M81 在 GALEX FUV 和 Spitzer-IRAC 4.5 hm 两 
个 波段 的 图 像 . 本 节 的 以 下 内 容 陈述 了 形态 学 解构 
前 的 初步 图 像 处 理 . 


2.2 ”银河 消光 改正 

首先 我 们 对 图 像 进 行 银 河 消光 改正 . 我 们 采用 
了 Fitzpatrick[ 约 在 1999 年 的 尘埃 消光 模型 , 其 中 选 
择 消光 与 总 消光 的 比 Rv 取 值 为 3.1，M81 在 V 波 段 
的 银河 消光 值 为 0.220 mag (Magnitude, 星 等 )[45]. 
银河 消光 改正 的 具体 过 程 由 名 为 extinction5 的 
Python 程序 包 完 成 . 


2.3 ”背景 噪声 扣除 

虽然 Spitzer 的 S4G 项 目 已 经 “官方 地 ”进行 了 
背景 噪声 扣除 , 但 是 他 们 扣除 的 背景 噪声 是 常数 
, 不 适用 于 背景 不 均匀 的 图 像 . 特别 是 , 面 亮度 径 
轮廓 拟 合 对 背景 非常 敏感 6, 因此 我 们 需要 构建 
一 个 二 维 背景 噪声 模型 用 于 背景 噪声 的 扣除 . 

背景 噪声 的 扣除 首先 需要 识别 背景 像素 . 我 
们 用 图 像 国 值 区 分 (Image Segmentation) 的 方法 ， 
以 3 倍 标准 偏差 、5 个 连接 像素 为 浆 值 来 识别 所 有 
信号 源 并 将 其 标记 出 来 , 这 样 就 对 图 像 中 信号 源 像 
素 和 背景 像素 做 了 区 分 . 

将 背景 像素 识别 出 来 以 后 , 所 有 背景 像素 都 将 
用 于 构建 二 维 背 景 模型 , 而 所 有 信号 源 像素 都 被 排 
除 在 背景 之 外 . 二 维 背景 模型 的 构建 利用 背景 像 
素 取 bin 的 方法 计算 局 部 变化 背景 的 尺度 大 小 中 1. 
该 方法 的 具体 操作 是 将 背景 像素 进行 取 bin, 通过 
计算 不 同 取 bin 窗 口 下 的 背景 像素 的 直方 图 来 确定 
局 部 背景 的 尺度 . 图 2 中 的 黑色 实 线 表示 原始 背景 
图 像 的 概率 分 布 直方 图 . 我 们 尝试 了 5-30 的 不 同 取 
bin O. 随 着 取 bin 窗 口 的 升 高 , 背景 像素 的 直方 图 


=} mt: 


[ 作 没有 选择 这 个 波段 
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将 会 变 得 越 来 越 狭窄 ; 同时 直方 图 的 均值 会 趋 于 稳 窗口 大 于 20 就 可 以 达到 直方 图 趋 于 稳定 的 效果 . 图 
E, 直到 不 再 变 狭 窄 , 因为 这 时 局 部 的 背景 像素 会 2 中 的 黑色 虚线 表示 取 bin 窗 口 为 20 的 背景 图 像 的 


) = 


全 部 被 取 bin 进 同一 个 像素 中 . 我 们 发 现 只 要 取 bin 概率 分 布 直方 图 . 


GALEX-FUV Spitzer-IRAC4.5um 
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图 1 M81 在 GALEX FUV ( 左 图 ) 和 Spitzer-IRAC 4.5 um ( 右 图 ) 的 图 像 ， 附 空间 坐标 以 展示 空间 尺度 . 


Fig.1 The GALEX FUV (left panel) and Spitzer-IRAC 4.5 um (right panel) images for M81, with coordinate superimposed for 


indicating the spatial scale. 
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图 2 M81 在 Spitzer-IRAC 4.5 pm 波段 的 原始 图 像 ( 左 图 )、 取 bin 后 的 图 像 (中 图 , 取 bin 窗 口 为 20)、 像 素 频率 分 布 直 方 统计 图 ( 右 图 ). 右 图 展示 


了 观测 图 像 的 背景 像素 (黑色 实 线 ) 与 取 bin 窗 口 等 于 20 的 背景 像素 (黑色 虚线 ) 的 像素 频率 分 布 直方 图 . 


Fig.2 The original (left panel) and binned (middle panel, with a 20-pixel window) images for M81 at the Spitzer-IRAC 4.5 um 
band. The right panel shows the pixel frequency distribution histogram of the background pixel in the observation image (black 
solid line) and the binned image with a binned window equal to 20 (black dotted line). 
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得 到 局 部 背景 的 尺度 大 小 以 后 , 我 们 将 观测 图 
像 分 割 成 多 个 矩形 的 子 区 域 , 这 些 子 区 域 的 大 小 
要 足够 包 进 前 景 星 , 但 是 也 要 小 于 局 部 背景 . 通过 
实验 我 们 取 相 对 保守 的 124 pixel x 104 pixel 的 区 
bk. 随后 将 每 个 这 样 的 区 域 合 并 成 一 个 像素 , 并 赋 
值 该 区 域 所 有 像素 的 3 倍 标准 偏差 的 中 值 , 其 中 我 
们 舍弃 掉 背 景 像素 数量 占 比 小 于 90% 的 区 域 , 随后 
形成 一 个 低 分 辩 率 的 二 维 背 景 图 像 . 最 后 用 双 三 次 


Original 


Background model 


图 3 M81 在 Spitzer-IRAC 4.5 hm 波段 的 观测 图 像 ( 左 


插值 , 将 该 低 分 辩 率 的 二 维 背 景 图 像 还 原 成 与 原始 
图 像 尺 寸 大 小 相同 的 背景 模型 图 像 , 随后 用 原始 图 
像 减 去 背景 模型 图 像 即 可 得 到 背景 扣除 后 的 数据 . 
在 此 我 们 将 背景 不 确定 性 定义 为 背景 扣除 产生 的 
不 确定 性 , 也 就 是 扣除 背景 后 的 所 有 背景 像素 值 的 
标准 偏差 . 背景 扣 减 结果 在 图 3 中 展示 . 本 工作 对 星 
系 M81 的 形态 学 解构 研究 基于 扣除 背景 后 的 图 像 
开展 . 


Backgound subtraction data 


辐 )、 背 景 模 型 图 像 (中 图 )、 扣 除 背 景 之 后 的 图 像 ( 右 图 ). 左 图 中 黑色 方 框 表示 用 于 背景 模 


型 建 模 所 


到 的 局 部 背景 的 区 域 , 黑色 区 域 是 没有 观测 的 区 域 . 为 了 清楚 显示 背景 的 空间 变化 , 我 们 增强 了 背景 图 像 的 对 比 度 . 


Fig.3 The observed (left panel), modeled background (middle panel), and background-subtracted (right panel) images for M81 


at the Spitzer-IRAC 4.5 um band. In the left panel, black boxes represent the area sampled for background modeling, and the 


black area was not observed. In order to clearly show the spatial variation in the background, we enhance the contrast scale for 


the background image. 


3 星系 图 像 形 态 学 解构 的 处 理 过 程 

在 本 工作 中 , 我 们 利用 GALFIT 软 件 P5, 28 gt 
行星 系 形态 学 解构 . 除了 需要 输入 待 解构 图 像 ， 
GALFIT 还 要 求 输入 屏蔽 图 像 (Mask Image). W 
Pi | (Sigma Image)j 以 及 点 扩散 函数 图 像 (Point 
Spread Function Image). 本 节 将 介绍 这 些 图 像 的 
制作 以 及 我 们 采用 的 星系 解构 模式 . 


屏蔽 、 噪 声 和 点 扩散 函数 图 像 的 制作 
屏蔽 图 像 用 来 将 某 些 像素 排除 使 之 不 参与 拟 
其 他 亮 源 以 及 图 像 中 的 坏 像素 对 星 
系 解构 产生 的 影响 . 被 屏蔽 的 目标 包括 前 景 星 、 非 
数字 信号 以 及 饱和 像素 等 . 


3.1 


前 景 星 是 主要 的 屏蔽 目标 . 我 们 的 工作 中 前 
景 星 分 为 两 类 : 一 类 是 M81 履 盖 天 区 之 外 的 前 景 星 
(不 与 星系 重 登 ), 另 一 类 是 与 M81 重 登 的 前 景 星 . 
此 探测 前 景 星之 前 需要 提前 确定 M81 履 盖 的 天 区 
范围 . 我 们 利用 等 照度 线 (Isophote) 来 确定 M81 禾 
盖 的 天 区 范围 . 在 拟 合 出 M81 的 椭圆 等 照度 线 之 
后 , 我 们 将 亮度 等 于 背景 不 确定 性 的 那个 标准 椭 
圆 等 照度 线 作 为 M81 的 边缘 , 该 椭圆 的 中 心 位 置 是 
Right Ascension (RA) = 148.888221?. Declination 
(Dec) = 69.065295*. ^E K fü 7507”. Æ G 加 
为 450”"、 方 位 角 为 -20”( 该 椭圆 也 是 后 续 对 星系 
进行 孔径 测 光 时 的 椭圆 孔径 ). 所 有 在 此 椭圆 以 儿 
的 前 景 星 将 会 用 图 像 闷 值 区 分 的 方法 来 标记 , 该 方 


LH 
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法 已 在 第 2.3 节 中 描述 . 的 噪声 图 像 还 需要 转换 成 与 观测 图 像 相同 的 单 
对 于 椭圆 孔径 内 的 前 景 星 , 我 们 采用 局 部 探测 us 

法 . 该 方法 类 似 于 第 2.3 节 中 采用 的 方法 , RIDEN PSF 图 像 在 GALFIT 对 星系 的 形态 学 解构 中 有 

到 内 部 分 成 多 个 子 区 域 ; 随后 我 们 在 每 个 子 区 域 设 。 “两 个 用 处 . 一 是 用 于 拟 合 图 像 中 的 点 源 , 二 是 对 模 


BS BUE, 将 高 于 该 闵 值 的 像素 定义 为 前 景 星 。 ”型 图 像 进行 卷 积 以 改变 其 空间 分 辨 率 . 


该 闵 值 根据 不 同 的 子 区 域 设置 ; 最 后 肉眼 检查 , 有 本 工作 中 PSF 图 像 制作 的 方法 是 尽 可 能 多 地 
必要 的 话 人 为 地 手动 添加 标记 和 去 标记 . 除了 前 景 找到 背景 扣除 后 图 像 中 的 没有 坏 数据 、 孤 立 (附近 
EMT STEREO DE Em. 没有 其 他 亮 源 污染 )、 相 对 不 太 亮 的 前 景 星 作为 样 


图 4 中 展示 了 去 除 所 有 前 景 星 的 M81 图 像 用 以 。。 本 , 通过 多 次 欠 代 制作 出 一 个 有 效 点 扩散 函数 (E 生 
检查 屏 项 效果 . 该 图 像 的 制作 方法 是 利用 观测 图 像 Da | hodie 
m mi ective Point Spread Function, ePSF) i79 Fk (& 49). 
TSR ER RD h EPDM RRL ee east J Ha ay BEAR eH, 


s 片 正 中 . ePSP 模 型 好 坏 可 以 通过 GALFIT 拟 合 来 检 


验 , 就 是 查看 GALFIT 是 否 能 只 用 一 个 ePSF 模 型 来 
很 好 地 拟 合 图 像 中 的 前 景 星 , 若 残 差 图 像 中 没有 该 
前 景 星 像素 的 明显 遗留 则 说 明 ePSF 模 型 可 靠 . 

图 5 中 展示 了 我 们 通过 以 上 方法 和 步 又 制作 上 
的 屏蔽 图 像 、 品 声 图 像 以 及 ePSF 图 像 . 具体 操 人 
由 名 为 photutils" 的 Python 软 件 包 实现 . 


3.2 ”M81 的 形态 学 解构 模式 


(1) Observation (2) Starclean 


TT LL 


图 4 M81 在 Spitzer-IRAC 4.5 hum 波段 的 观测 图 像 ( 左 图 ) 和 去 除 前 本 节 的 内 容 包 括 描述 在 星系 解构 过 程 中 需要 
景 星 后 的 图 像 ( 右 图 ) 的 展示 . 两 张 图 的 空间 尺度 相同 . ç : " FREUEN "m 
` 拟 合 的 子 结构 、 如 何 选择 函数 、 如 何 建 模 以 及 模 

Fig.4 Observed (left panel) and star-cleaned images (right J AENT EN 4 E Fa Et ae 站 的 其 

panel) for M81 at the Spitzer-IRAC 4.5 um band. The two 型 好 坏 的 Vt. 核 球 和 a T AUS J 的 基本 

images have the same physical scale. 子 结 KJ, 他 们 在 星系 中 的 质 量 和 光度 的 占 比 非 常 
高 , 很 大 程度 上 决定 了 星系 的 形成 和 演化 . 本 工作 


GALFIT 的 拟 合 原理 是 最 小 一 乘法 , ETATE 。 ”的 科学 目标 是 探究 M81 的 子 结构 , 并 且 对 这 些 子 结 
声 图 像 参与 最 小 二 乘 计算 . 噪声 图 像 包含 泊 松 噪声 构 的 形态 参数 进行 量化 . 工作 中 我 们 把 核 球 和 盘 作 


和 读 出 噪声 , 遵循 以 下 公式 : 为 M81 的 基本 子 结构 , 在 不 同 的 解构 模式 中 , 我 们 
添加 不 同 的 函数 来 拟 合 其 他 子 结构 . 多 种 解构 模式 


o° = op + ORo (4) 
有 助 于 调查 不 同 子 结构 的 拟 合 对 核 球 和 盘 的 形态 
= V+. G 学 参数 结果 的 影响 
其 中 表示 输入 的 噪声 图 像 , op 表示 泊 松 噪声 , on GALFIT 中 内 置 了 许多 函数 , 包括 Sérsic 函 数 、 


表示 读 出 噪声 , 到 表示 以 电子 为 单位 的 未 扣除 背景 。 ”指数 函 数 (指数 函数 在 GALFIT 中 简称 为 “Expdisk” 
的 观测 图 像 6. 计算 泊 松 噪声 的 以 电子 为 单位 未 扣 函数 , 下 文 我 们 用 'EExpdisk” 表 示 对 盘 结 构 进 行 拟 
除 背 景 的 图 像 需 要 通过 图 像 文件 头 (Header) 中 的 合 的 指数 函数 )、 高 斯 函数 、 点 扩散 函数 等 , 这 些 
GAIN (CCD 的 增益 ) 参 数 来 得 到 . 读 出 噪声 函数 被 称 为 “基础 模型 函数 ", 其 中 S6rsic 函 数 是 最 
由 Header 中 的 RONOISE 参 数 记录 . 由 此 公式 计算 出 常用 的 描述 系统 径 向 面 亮度 轮廓 的 函数 , 函数 形式 


Shttps://users.obs.carnegiescience.edu/peng/work/GALFIT/GALFIT.html. 
"https:/ /photutils.readthedocs.io/en/stable/index.html. 
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如 (3) 式 , 由 于 其 较 高 的 自由 度 因此 可 以 广泛 地 用 于 
不 同 子 结构 包括 但 不 限于 核 球 、 盘 、 星 系 棒 等 的 建 
模 ; Expdisk 函 数 是 Skrsic 函 数 的 SErsic 指 数 等 于 1 的 
特殊 形式 ; 高 斯 函数 是 Srsic 函 数 的 Skrsic 指 数 等 于 
0.5 的 特殊 形式 . 基础 模型 函数 用 来 拟 合 面 亮 度 均 
匀 规 则 的 子 结构 . 为 了 拟 合 旋 臂 、 环 等 面 亮度 不 
均匀 、 不 规则 的 子 结构 , GALFIT 内 置 了 分 别 为 方位 


(1) Observation 


函数 (Azimuthal Function) 4r RK Zi (Truncation 
Functions) 的 两 种 “高 级 函数 ”. 这 两 种 高 级 函数 不 
能 单独 使 用 , v BE Bt DON ES a 其 
作用 是 将 基础 函数 修改 成 特殊 的 不 规则 形状 . 方位 
函数 中 包含 一 个 坐标 旋转 函数 , 通过 将 其 登 加 在 基 
础 函数 上 可 以 产生 旋 臂 的 结构 ; 截断 函数 通过 截断 
基础 函数 从 而 产生 环 结构 或 者 断裂 的 效果 . 


(3) Sigma image (4) PSF image 


* 


图 5 M81 在 Spitzer-IRAC 4.5 hm 波段 的 观测 图 像 (1)、 屏 项 图 像 (2)、 噪 声 图 像 (3)、PSE 图 像 (4), 其 中 图 (1)、(2)、(3) 的 物理 尺度 相同 . 


Fig.5 Observed (1), mask (2), sigma (3), and PSF (4) images for M81 at the Spitzer-IRAC 4.5 um band, where the physical 


scales in the images (1), (2), and (3) are the same. 


M8108 T I PR FU EX WW #h f 25 J ULP, 还 有 
两 条 非常 宏伟 壮观 的 旋 辟 ; 星系 中 心 也 有 一 个 
AGN[59,; 此 外 , 核 球 中 存在 一 个 若隐若现 的 类 似 旋 
辟 的 子 结构 [51. 我 们 将 会 对 这 些 子 结构 依次 建 模 . 
我 们 选择 自由 度 高 也 是 最 常用 的 S6rsic 函 数 作为 以 
上 所 有 子 结构 的 基础 函数 , 另外 在 Srsic 函 数 的 基 
础 上 全 加 坐标 旋转 函数 或 者 截断 函数 来 构建 旋 辟 
结构 、 用 PSF 函 数 来 拟 合 星 系 中 心 AGN 结 构 . 有 很 
多 工作 都 是 利用 Expdisk 函 数 来 拟 合 星系 盘 , 我 们 
在 解构 中 没有 遵循 这 个 传统 , 因为 盘 的 SGrsic 指 数 
也 是 我 们 想 探 寻 的 参数 , 特别 是 检查 不 同 子 结 构 拟 
合 对 盘 的 集中 度 的 影响 ， 

所 有 这 些 模型 在 图 像 中 的 坐标 遵循 椭圆 的 分 
布 . 在 拟 合 旋 臂 时 , 坐标 旋转 函数 用 于 将 原本 按照 
椭圆 分 布 的 每 个 像素 的 坐标 强制 按照 新 的 方式 重 
新 分 布 , 相当 于 对 原本 的 椭圆 函数 进行 扭曲 , 因此 
可 用 来 拟 合 不 同形 态 的 旋 臂 结构 . M81I 的 外 旋 臂 
在 IRAC 4.5 hm 波段 的 亮度 相对 比 于 核 球 和 盘 来 
说 并 不 太 高 . 为 了 保证 构建 出 合理 的 外 旋 臂 模型 ， 


我 们 利用 旋 臂 结构 异常 明显 的 GALEX-FUV 图 像 
作为 参照 . 从 图 1 中 可 以 看 到 外 旋 臂 并 不 与 核 球 直 
接 接壤 , 而 是 旋 辟 内 部 存在 明显 的 “ 空 区 ”, 因此 我 
们 在 构建 外 旋 臂 的 时 候 需 要 考虑 截断 . 在 拟 合 M81 
的 外 旋 臂 时 , 我 们 将 坐标 旋转 函数 和 截断 函数 分 
别 县 加 在 Sérsic 函 数 上 来 构建 旋 辟 结构 . 我 们 先 在 
GALEX-FUV 上 拟 合 外 旋 恬 模型, 然后 把 这 个 模型 
应 用 在 IRAC 4.5 pm 图 像 进行 拟 合 . 
本 工作 中 的 坐标 旋转 函数 的 形式 为 究 律 双 曲 
正切 坐标 旋转 , 函数 表达 式 如 下 : 


6(r) = Oout tanh (ris; Tout; bined; gu r) x 


1 T s 
EM (6) 
E rp, O(r) #2 E fer Mb AY e pz f BE; mm 表示 旋 臂 


开始 的 径 向 距离 , 而 rou 表 示 旋 臂 结束 的 径 向 距离 ; 
当 r 小 于 rouw 时 旋 臂 形式 主要 由 tanh 函 数控 制 的 , Fe 
律 相 控制 的 是 rou 之 外 旋 臂 特征 ,oa 是 桶 律 斜 率 ; 
bu 表示 旋 臂 结束 的 径 向 距离 rou 处 的 旋转 角度 ， 
因此 旋 臂 的 缠绕 松紧 程度 由 gs 和 rout 共 同 控制 ; 
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zr Te E ana BOE EA rin EO (m) LP A E A, 


常用 于 拟 合 星系 棒 ; 对 于 无 棒 星 系 , ri 常常 被 设置 
为 0; 旋 臂 盘面 的 倾斜 角度 是 ga, 而 ba 是 方位 角 ， 


这 两 个 参数 共同 定义 模型 在 天 球 上 的 投影 . 
本 工作 中 的 截断 函数 的 数学 表达 式 为 : 


P(x, y) = tanh(z, y; zo; Yo, Tbreak, AT softs q, OP.A.) ; 

(7) 
其 中 P(x,y) 表 示 在 坐标 (zx,y) 处 的 截断 因子 , 截断 
函数 通过 分 别 将 截断 因子 PP 或 1 — PP 与 截断 内 部 或 
戌 断 外 部 的 原始 流量 分 布 相 乘 来 截断 模型 . (o, yo) 
是 中 心 坐 标 , g 是 轴 比 , Op. A 是 截断 函数 的 方位 角 
(PA); rpreark 是 断裂 半径 , 表示 截断 模型 流量 到 该 半 
径 处 会 下 降 至 其 原始 流量 99%; Arsog 为 缓冲 长 度 ， 
被 定义 为 rsof — rbreak 或 rureak — som TAIA FARR 


SUIT EX, P EIUS, 其 中 rsok 是 被 截断 的 模型 流量 
降 到 其 原始 流量 的 1% 的 半径 . 

除了 以 上 用 于 拟 合 延展 的 结构 的 函数 以 外 ， 
PSF 函 数 是 最 常用 于 拟 合 点 状 结构 的 函数 . 不 过 
不 同 数 据 可 能 表现 出 完全 不 同 的 PSF 形 状 , 因此 
GALFIT 需 要 用 户 自行 提供 PSF 图 像 来 拟 合 点 源 结 
构 或 者 用 于 图 像 的 卷 积 . 本 工作 中 的 PSF 的 模型 图 
像 来 自在 第 3.1 节 中 制作 的 ePSF 图 像 . 

在 开始 星系 解构 之 前 , 需要 人 为 提供 一 个 模 
WME, GALFIT 软 件 需 要 读 取 该 文件 后 才能 运行 . 
模板 文件 包括 两 部 分 , 分 别 是 限制 拟 合 过 程 的 参 
数 和 模型 函数 的 初始 猜测 值 . 对 于 GALFIT 需 要 人 


为 地 输入 一 些 限制 拟 合 过 程 的 参数 , 例如 拟 合 窗 
口 大 小 、PSF 卷 积 盒 大 小 等 , 所 有 的 这 些 参数 我 


们 都 会 根据 GALFIT 官 网 的 Frequently Asked tech- 
nical Questions (FAQ) 中 的 建议 来 设置 . 其 中 拟 合 
窗口 和 PSF 卷 积 核 的 大 小 很 大 程度 上 可 能 影响 拟 
合 时 间 和 拟 合 的 准确 性 , 原则 上 拟 合 窗口 和 PSF 卷 
只 核 窗 口 大 小 都 要 尽 可 能 大 , 以 便 完全 容纳 目标 
星系 和 PSF 的 所 有 像素 ; 但 是 太 大 的 拟 合 和 卷 积 核 
窗口 会 降低 拟 合 的 速度 , 因此 选择 合适 的 拟 合 和 
卷 积 核 窗 口 能 提高 拟 合 的 效率 . 我 们 的 图 像 有 足 
够 大 的 视 场 , 因此 拟 合 窗口 我 们 选择 2000 pixel x 
2000 pixel, 这 个 尺度 足以 包括 M81 全 部 以 及 很 大 
一 部 分 背景 区 域 ; 卷 积 核 设置 80 倍 的 图 像 PSF 的 半 


20-8 


期 


高 全 宽 , 像素 尺度 和 测 光 零点 遵循 Spitzer 望 远 镜 手 
册 设 置 为 0.75” /pixel 和 18.32 mag. 

在 开始 进行 形态 学 解构 拟 合 之 前 , 还 需要 设 定 
模型 的 初始 参数 . 初始 参数 对 于 相对 简单 、 没 有 用 
到 坐标 旋转 和 截断 函数 的 模型 来 说 影响 不 大 , 特别 
是 对 于 单 Sérsic 函 数 的 拟 合 , 即使 初始 参数 的 设置 
极端 不 合理 , 最 后 都 能 稳定 地 得 到 同一 个 结果 . 
而 , 对 于 坐标 旋转 和 截断 函数 , 初始 参数 的 设置 会 
对 GALFIT 的 拟 合 产 生 影响 . 以 我 们 在 本 工作 中 得 到 
的 经 验 来 说 , 初始 参数 影响 坐标 旋转 和 截断 函数 拟 
合 的 时 间 长 短 , 合理 的 参数 猜测 通常 带 来 更 短 的 拟 
合 时 间 . 对 于 参数 的 设置 , 我 们 采取 一 个 相对 经 验 
的 方法 来 估计 所 有 函数 的 初始 参数 值 : 在 输入 初始 
参数 前 人 为 创建 一 个 和 目标 看 起 来 相像 的 人 工 模 
型 . 这 个 人 工 模 型 的 构建 是 通过 人 为 手动 调节 参数 
得 到 , 通过 输出 该 图 像 和 视觉 检查 , 确保 人 工 模型 
图 像 和 目标 相似 . 例如 人 工 模型 比 目 标 更 亮 就 降低 
亮度 , 人 工 模型 比 目 标 小 就 升 高 有 效 半 径 等 . 接着 
尝试 用 这 个 人 工 模型 的 参数 值 作为 初始 参数 的 设 
EL. 由 于 这 些 复杂 的 不 规则 模型 的 先 验 条 件 对 拟 合 
有 更 大 的 影响 , 因此 人 工 输出 一 个 好 模型 的 目的 是 
提供 一 个 好 的 先 验 条 件 . 需要 注意 的 是 , 除了 外 旋 
臂 是 根据 GALEX-FUV 图 像 进 行人 工 模型 建 模 的 
以 外 , 其 他 子 结构 人 工 模型 的 建 模 都 是 在 Spitzer- 
TRAC 4.5 hm 的 图 像 上 进行 . 由 于 旋 臂 模型 需要 在 
Sérsic 函 数 的 基础 上 县 加 坐标 旋转 或 者 截断 函数 ， 
自由 参数 要 比 其 他 能 通过 均匀 单 Strsic 函 数 描述 的 
模型 更 多 , 因此 旋 臂 模型 的 参数 设置 需要 花费 比 核 
球 和 盘 这 些 规则 模型 更 多 的 时 间 和 精力 , 需 反 复 实 
验 以 保证 得 到 一 个 足够 真实 的 拟 合 结果 . 为 了 保证 
所 有 解构 策略 的 相对 独立 性 , 不 同 策略 的 函数 初始 
参数 设置 都 同样 按照 上 述 方法 来 估计 , 也 就 是 说 我 
们 会 假设 每 次 进行 新 的 星系 解构 策略 的 实验 之 前 
都 会 假设 这 是 第 1 次 实验 . 

本 工作 中 , 我 们 采用 经 验方 法 来 估计 GALFIT 的 
解构 拟 合 误差 (6, 92, 我 们 在 待 解 构图 像 上 分 别人 
工 加 上 或 减 去 1 倍 背 景 不 确定 性 的 常数 值 后 重新 拟 
合 , 新 的 拟 合 结果 将 会 作为 误差 的 上 限 和 下 限 , 这 
个 方法 应 用 于 除了 星 等 以 外 的 所 有 参数 的 误差 估 
th. 星 等 (或 亮度 ) 估 计 的 不 确定 性 来 源 不 只 是 背景 
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的 不 确定 性 , WR FL 128 v ZR ER pu B DC IJ E VA 
差 都 是 亮度 不 确定 性 的 来 源 . 因此 我 们 将 椭圆 孔 径 
的 误差 、 定 标 误差 和 背景 不 确定 性 的 正 交 组 合 ( 即 
平方 和 的 开 方 ) 作 为 星 等 (或 亮度 ) 的 误差 . 


结果 展示 
本 工作 中 , 我 们 采用 由 简单 到 复杂 的 6 种 模式 
( 即 6 种 子 结构 组 合 ) 在 Spitzer-IRAC 4.5 pm 图 像 上 
对 星系 进行 解构 : 


e 单 成 分 : Sérsic (M81 整 体 ) 
e 双 成 分 : Sérsic (M81 整 体 )+PSF (AGN) 
e IURI 


e 三 成 分 
(AGN) 


e 四 成 分 : Sérsic (KER) + Sérsic (4) + PSF 
(AGN) + Coordinate rotation-Truncated- 
Sérsic (外 旋 臂 ) 


e ERI: Sérsic ( 核 球 ) + Sérsic (4) + PSF 
(AGN) + Coordinate rotation-Truncated- 
Sérsic (外 旋 臂 ) + Coordinate rotation-Sérsic 
(内 旋 臂 ) 


对 于 每 种 模式 , 我 们 都 会 实际 测量 星系 整体 以 
及 各 子 结构 成 分 的 面 亮度 径 向 轮廓 , 用 来 和 模型 进 
行 比较 . 面 亮 度 径 向 轮廓 的 测量 方法 是 对 已 去 除 前 
景 星 的 M81 图 像 ( 见 图 4), 画 出 等 照度 线 来 进行 对 面 
亮度 的 采样 , 所 有 等 照度 线 的 椭圆 与 M81 测 光 孔 径 
的 中 心 、 轴 比 和 位 置 角 相 等 . 这 组 等 照度 线 也 统一 
用 于 所 有 子 结构 图 像 . 图 6-11 分 别 展示 了 每 种 解构 
模式 的 结果 图 像 和 面 亮度 径 向 轮廓 , 图 12 是 各 模型 
TASSE 合 结果 展示 . Aper 吉 构 的 形 
态 学 参数 值 记录 在 表 1 中 . 以 下 各 小 节 分 别 旦 现 每 
种 模式 的 结果 分 析 . 


41 单 成 分 拟 合 (IM81 整 体 ) 


只 用 单 S6rsic 函 数 进行 的 单 成 分 拟 合 得 到 M81 
的 整体 Sérsic 指 数 为 4.29. 原始 图 像 和 模型 的 测 光 


4 


— 


: Sérsic ( 核 球 ) + Sérsic ($) 


: Sérsic (FE ER) + Sérsic ($) + PSF 
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结果 显示 , 单 成 分 拟 合 的 模型 星 等 与 在 无 前 景 
星 图 像 上 进行 孔径 测 光 得 到 的 结果 一 致 (6.78 + 
0.02 mag). 表 1 中 的 光度 占 比 的 定义 为 在 相同 孔径 
下 子 结构 模型 图 像 或 者 所 有 子 结构 模型 总 和 的 模 
型 图 像 在 原 星系 总 流量 中 的 占 比 , 由 于 模型 可 能 

不 同 程 度 地 高 估 或 低估 星系 的 亮度 , 因此 所 有 子 
结构 的 占 比 总 和 不 严格 等 于 100%. 在 单 成 分 拟 合 
的 结果 中 , M81 整 体 的 模型 略微 高 佑 了 M81 的 亮度 ， 


因此 其 光度 占 比 为 100.51%. 残 差 图 像 的 光度 比 为 
12.7396, 主要 由 残留 的 旋 臂 结构 、 少 量 的 核 球 、 盘 
以 及 其 他 恒星 提供 . 

然而 , 如 图 6 所 示 , 只 用 单 S6rsic 函 数 来 拟 合 是 
一 种 粗略 用 于 测量 星系 整体 性 质 的 方法 . 从 图 中 可 


以 看 出 , 模型 图 像 低 估 了 星系 的 方位 角 ; 残 差 图 像 
显示 模型 高 估 了 星系 在 小 半径 区 域 以 及 旋 臂 间 的 
亮度 , 低估 了 旋 臂 以 及 星系 盘 上 的 亮度 , 因此 导致 
了 旋 臂 间 和 部 分 旋 臂 被 过 度 拟 合 而 残 差 图 像 还 剩 
余 了 不 少 的 子 结构 . 面 亮 度 径 向 轮廓 图 清晰 地 显 
示 单 成 分 拟 合 模 型 的 椭 率 和 方位 角 只 是 近似 等 于 
观测 值 的 平 值 , 而 无 法 准确 地 反映 椭 率 和 方位 角 
随 半 径 的 变化 趋势 . 此 外 , 模型 明显 高 佑 了 M81 在 
400” 以 外 的 亮度 , 并且 随 着 半径 的 增加 , 模型 和 
观测 图 像 的 匹配 程度 也 逐渐 下 降 . 这 也 说 明 单 成 
分 Skrsic 函 数 在 估计 星系 的 面 亮 度 轮廓 时 有 缺陷 ， 
因为 星系 在 形态 学 上 的 复杂 性 导致 星系 整体 的 面 
亮度 径 向 轮廓 并 非 严 格 遵循 某 个 特定 函数 的 形式 . 
单 成 分 拟 合 M81 的 形态 学 成 分 是 一 个 过 于 简单 的 
尝试 . 

4.2 ” 双 成 分 拟 合 (M81 整 体 十 AGN) 

M81 中 心 存在 一 个 低 光度 的 AGN[59. 与 寄主 
星系 的 辐射 来 源 不 同 , AGN 的 辐射 不 属于 任何 
一 种 星 族 ; 但 是 AGN 的 存在 可 能 会 对 星系 的 动力 
学 、 形 态 学 、 物 理性 质 甚至 其 形成 和 演化 产生 影 
响 53-55]; 就 观测 而 言 , 如 果 不 扣除 AGN 的 贡献 , 那 
么 寄主 星系 的 观测 特征 也 会 因 AGN 的 存在 而 产生 
偏差 . 因此 , 如 何在 测量 星系 的 时 候 扣 除 AGN 的 贡 


献 一 直 都 是 星系 研究 中 一 个 重要 的 问题 . 本 工作 
中 , 我 们 尝试 利用 星系 形态 学 解构 将 AGN 和 寄主 
星系 分 离 . 
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3948326. 早 型 旋涡 星系 M81 (NGC 3031) 的 形态 学 解构 研究 2 期 


图 6 单 成 分 的 拟 合 结果 展示 . (1) 表 示 原 始 图 像 ; (2) 表 示 模 型 


包括 模型 和 原始 图 像 的 椭 率 (1 一 b/a, 其 


(4) Surface brightness profile 


(3) Residual 


1 


01 


107 
Radius/” 


Pa. bU AS ERE 9 
度 轮廓 中 的 圆圈 表示 观测 数据 ， 黑 


图 像 ; (3) 表 示 残 差 图 像 ， 由 观测 图 像 碱 去 模型 图 像 获得 ; (4 表示 径 向 面 亮度 轮廓 图 
长 轴 和 


FA 


轴 ), 方位 角 、 面 亮度 (以 mag/(”)” 为 单位 ) 和 面 亮度 残 差 的 比较 . 面 亮 
色 虚 线 表示 模型 图 像 的 数据 . 


Fig.6 Results of single-component fitting. The panel (1) represents the observed image; the panel (2) represents the model 


image; the panel (3) represents the residual image, which is obtained by subtracting the model image from the observed image; 


the panel (4) represents the radial surface brightness profile, including the comparison of ellipticity (1 — b/a, where a and b are 


semimajor axis and semiminor axis of ellipse respectively), the position angle, surface brightness (in unit of mag/(")?) and 


surface brightness residual between the model and the observed images. The open circles in the surface brightness profiles 


represent the observational data, and the black dotted lines represent the model curves. 


RAJH Sérsic K APSF A Z H 
果 如 图 7 所 示 , 子 结构 图 像 在 图 12 


体 以 及 AGN, 其 结 


中 , 子 结构 参数 值 记 录 在 表 1 中 . 图 7 显示 模型 
和 残 差 图 像 与 单 成 分 拟 合 结果 相差 不 大 ; 径 向 
度 轮廓 也 与 单 成 分 拟 合 的 结果 相近 , 都 是 模型 
无 法 精准 地 描述 椭 率 和 方位 角 随 半径 的 变化 趋势 ; 


| 拟 合 M81 整 


图 像 


同样 , BER AY h 


分 的 亮度 . 因为 M81 的 AGN 亮 度 低 , 所 以 对 寄主 星 


系 的 影响 很 弱 . 


亮度 高 估 了 距 星 系 中 心 400” 以 


图 像 


外 部 


在 表 1 中 , AGN 比 M81 整 体 暗 弱 很 多 ( 星 等 差 为 
6 mag). 在 添加 了 AGN 成 分 之 后 M81 整 体 的 有 效 半 


径 上 升 了 24.3%%; Sérsict Be F £ Y 9.32%; 轴 


比 和 
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方位 角 几 乎 无 变化 . 添加 了 AGN 模 型 之 后 , M81% 
体 的 光度 占 比 下 降 了 0.79%, 而 AGN 本 身 的 亮度 只 
占 总 亮度 的 0.64%. 变化 最 大 的 是 有 效 半径 , 这 是 因 
为 考虑 了 AGN 模 型 之 后 M81 整 体 模型 亮度 的 重新 
分 配 导 致 的 , 但 是 变化 范围 比 背 景 不 确定 性 带 来 的 
误差 要 更 小 . 残 差 图 像 的 光度 比 为 12.38%, 与 单 成 
分 拟 合 结果 相似 . 


4.3 ” 双 成 分 拟 合 ( 核 球 十 盘 ) 

对 于 同时 具有 核 球 和 盘 结 构 的 星系 , 人 们 常用 
核 球 上 盘 的 双 成 分 来 进行 形态 学 解构 56-584. M81 就 
是 这 样 的 星系 , 因此 核 球 + 盘 的 双 成 分 形态 学 解构 
(两 个 Sérsic 函 数 ) 也 被 我 们 用 于 6 种 模式 之 一 . 
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(3) Resid 


Ka 


7 双 成 分 (M81 整 体 +AGN) 的 拟 合 结果 持 


ual 


Radius/” 


ZN. 


到 像 的 格式 和 


图 6 相同 . 面 亮度 轮廓 中 不 同 的 线 表 示 子 结构 模型 的 数据 . 


Fig.7 The fitting results of two components (M81+AGN). The format of the image is the same as that in Fig. 6. Different lines 


in the surface brightness profile represent the data of the substructure model. 


图 8 展示 了 双 成 分 ( 核 球 + 盘 ) 解 构 的 图 像 结果 ， 
其 中 各 个 成 分 的 模型 图 像 结果 可 见 图 12. 相 比 单 


Sérsic 函 数 的 拟 合 结果 , 双 成 分 ( 核 球 + 盘 ) 拟 合 能 够 


更 好 地 描述 M81 整 体 和 局 部 的 形态 . 从 模型 图 像 的 
大 小 、 轴 比 和 方位 角 都 非常 接近 观测 图 像 . 残 差 图 


像 显示 模型 很 好 地 将 核 球 和 盘 从 观测 


图 像 中 抽 离 


出 来 了 , 只 


剩 下 了 两 条 明亮 外 旋 臂 ; 星系 内 部 的 内 
旋 臂 以 及 星系 中 心 类 似 小 型 星系 棒 的 结构 残留 在 
图 像 上 . 同时 模型 的 椭 率 和 方位 角 的 径 向 变化 
趋势 与 观测 结果 比较 一 致 . 相 比 此 前 单 S6rsic 函 数 


在 半径 超过 400” 以 后 径 向 面 


描述 M81 形 态 的 原因 


亮度 轮廓 估 值 偏 高 的 


情况 也 有 所 改善 . 核 球 和 盘 的 大 小 、 椭 率 和 方位 角 
都 有 明显 的 不 同 , 这 也 是 单 S6rsic 函 数 不 


能 很 好 地 


ERIP, 核 球 和 盘 之 间 
大 约 0.3 星 等 的 差 值 ; 从 亮度 


的 亮度 差异 不 大 ， 只 J 
占 比 来 看 , 核 球 和 盘 分 


别 占据 了 42.00% 和 56.81%% 的 流量 , 两 个 成 分 亮度 差 


别 不 是 很 大 . 这 说 明 即 使 是 如 


E 主 要 探测 年 老 恒 星 辐 


射 的 4.5 pm 的 波段 上 , 稳 的 亮度 也 是 不 可 忽视 的 . 
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核 球 的 有 效 半 径 (re) 仅 为 盘 的 约 1/3, 核 球 的 Sérsic 


BB(n = 5.18) 要 远大 于 盘 
说 明 昌 然 核 球 和 盘 的 亮度 相近 , 但 是 核 球 集中 在 星 


的 Sérsic 指 数 (n = 1.17), 


系 小 半径 区 域 中 , 而 盘 的 分 布 更 为 广阔 和 平坦 . 核 


球 和 盘 的 轴 比 和 方位 角 有 显著 的 不 同 , 这 也 是 径 向 


Ti > BE se Eg E| rH 3 
之 一 . W^ Sérsic] 


比 和 方位 角 随 半 径 变 化 的 原因 
数 拟 合 的 模型 


图 像 同样 可 以 在 


误差 范围 内 精准 地 预测 M81 的 亮度 . 我 们 发 现 残 差 


图 
内 旋 臂 和 其 


4.4 


果 基 本 一 致 ; AGN 只 
弱 很 多 , 其 光 


度 占 比 约 为 0.13%. 残 差 图 


像 的 光度 占 比 为 10.12%, 主要 由 剩余 的 外 旋 臂 、 
他 恒星 提供 . 


三 成 分 拟 合 ( 核 球 十 盘 十 AGN) 
接着 我 们 在 核 球 + 盘 的 双 成 分 基础 上 添加 
AGN 成 分 , 得 到 的 结果 见 图 9 和 12 以 及 表 1. 这 种 
模式 下 核 球 和 盘 的 结果 与 单纯 两 个 Srsic 函 数 的 结 
TRAMES, 比 核 球 和 盘 暗 


像 与 两 个 


SeErsic 函 数 的 结果 相 比 也 没有 明显 变化 , 其 光度 占 


比 约 为 10.10%. 
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(4) Surface brightness profile 


SESS Sacs * Bulge--Sérsic 
?SSseo、、 = = Disk--Sérsic 


(3) Residual 


Radius/” 


8 双 成 分 ( 核 球 十 盘 ) 的 拟 合 结果 展示 . 图 像 的 格式 和 图 6 与 图 7 相同 . 


Fig.8 The fitting results of two components (bulge + disk). The format of the image is the same as that in Fig. 6 and Fig. 7. 


(1) Observation 
ora (4) Surface brightness profile 


= Model 

= * Bulge--Sérsic 

sssss= = = Disk--Sérsic 
sSescoo Se O Data 

*  AGN-PSF 


(3) Residual * 
< e he e Ae He e Ae He He He e de He de He de He He de de ke he kehe ke keke kke k kek kk kek 
$ maaana a a oo Se JUD. 


Radius/" 


图 9 三 成 分 的 拟 合 结果 展示 . 图 像 的 格式 和 图 6 与 图 7 相同 . 


Fig.9 The fitting results of three components. The format of the image is the same as that in Fig. 6 and Fig.7. 
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4.5 ”四 成 分 拟 合 ( 核 球 十 盘 十 AGN 十 外 旋 辟 ) 

之 前 的 结果 都 显示 在 残 差 图 像 中 有 明显 的 旋 
臂 结构 , 因此 在 上 一 节 解 构 模 式 的 基础 上 增加 了 外 
旋 辟 成 分 , 即 核 球 + 盘 +AGN+ 外 旋 臂 的 四 成 分 解 
构 模 式 ( 从 之 前 的 结果 中 看 , M81 的 旋 臂 似乎 有 不 
连续 的 内 外 两 段 构成 , 这 一 节 的 解构 我 们 只 考虑 外 
旋 臂 , 下 一 节 我 们 将 加 入 内 旋 臂 )， 

与 其 他 子 结构 不 同 , Je E i BE 
指数 上 县 加 坐标 旋转 函数 (6) 式 和 截断 函 
因此 旋 臂 的 自 
初始 参数 的 设置 对 结果 的 影响 更 大 ; 
近 红 外 波段 , 旋 臂 可 能 会 被 其 他 子 结 
导致 初始 参数 的 设置 更 加 困难 . 因 


| 
[n 


N 


EJ3 5] Sérsic 
数 (7) 式 ， 
由 参数 比 均匀 平滑 的 结构 更 多 ， 
司 时 , 在 

KJ AE YZ, 
此 我 们 利 


jGALEX-FUV 的 图 像 来 确定 旋 辟 的 初始 参数 . 


具体 做 法 是 ， 人 工 调 节 各 参数 产生 一 个 与 FUV 图 
像 一 致 的 模型 ， 再 用 这 些 参 数值 通过 简单 的 


单位 换算 作为 在 Spitzer-IRAC 4.5 hm 图 像 拟 合 
时 的 初始 参数 (其 中 旋 臂 的 亮度 设置 为 比 核 
球 高 2 个 星 等 ). 因为 外 旋 臂 距离 星系 中 心 大 约 
100” 处 截断 , 所 以 在 外 旋 辟 的 建 模 时 需要 加 入 截 
Wr eR C 


核 球 + 盘 +AGN+ 外 旋 臂 的 四 成 分 解构 模式 的 


结果 见 图 10 和 12 以 及 表 1. 从 图 10 和 12 中 可 以 看 出 ， 
在 添加 了 外 旋 辟 之 后 , 模型 图 像 与 观测 图 像 的 相似 
度 有 明显 提高 ; 残 差 图 像 中 也 几乎 看 不 到 外 旋 臂 的 
Pm, 图 像 中 除了 一 些小 尺度 的 亮 恒星 形成 区 以 外 
Ul EE. 模型 的 面 亮 度 径 向 轮 廊 也 能 更 好 地 拟 
合 观测 ; 模型 的 椭 率 和 方位 角 随 半径 变化 的 曲线 
也 与 观测 更 加 一 致 , 同时 , 此 前 的 模型 对 星系 外 围 
面 亮 度 的 高 估 也 有 非常 大 的 改善 , 甚至 到 了 接近 
700” 的 位 置 , 模型 与 观测 图 像 的 面 亮 度 差 值 也 只 有 
约 0.6 星 等 . 对 比 单个 子 结构 , 添加 外 旋 辟 模型 前 后 
对 核 球 模型 图 像 没 有 明显 的 影响 ; 但 是 添加 了 外 旋 
辟 模 型 之 后 , 盘 的 大 小 发 生 了 轻微 的 变化 . 

从 表 1 中 可 以 看 出 , 旋 辟 的 添加 没有 对 核 球 和 
盘 产 生 非常 大 的 影响 , 大 多 数 参数 受到 的 影响 比 背 
景 误差 更 小 . 变化 最 大 的 参数 是 盘 的 有 效 半径 (re)， 
添加 外 旋 辟 后 盘 的 x。 下 降 约 11%. 此 外 , 外 旋 辟 的 添 


= 
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加 会 使 核 球 和 盘 的 方位 角 发 生 细 微 的 变化 . 从 光度 
看 , 外 旋 辟 的 光度 只 占星 系 整体 的 5.42%, 添加 外 旋 
臂 后 的 核 球 和 盘 的 光度 比 都 有 略微 下 降 , 而 AGN 没 
有 太 大 变化 . 残 差 图 像 的 光度 比 约 为 9.15%, 主要 来 
自 内 旋 臂 和 其 他 恒星 


4.6 ”五 成 分 拟 合 ( 核 球 十 盘 十 AGN 十 外 六 

WEE) 

在 以 上 的 结果 中 都 能 看 到 M81 存 在 一 个 内 旋 
B. 很 多 对 M81 的 研究 都 发 现 了 这 个 与 外 旋 臂 不 连 
续 的 子 结构 6, 597601, 这 些 研究 结果 表明 , M81 的 内 
旋 辟 和 外 旋 避 都 遵循 这 同一 个 旋转 模式 , 但 是 组 成 
物质 可 能 不 同 . 因为 内 旋 臂 处 在 小 半径 区 域 , 与 核 
球 直接 接壤 , 所 以 如 果 不 考虑 该 子 结构 , 那么 在 解 
构 中 来 自 内 旋 臂 的 光 可 能 会 被 分 配 到 核 球 或 者 盘 
E, 导致 在 结果 中 产生 误差 . 因此 在 这 一 节 中 我 们 
把 这 个 结构 添加 进 拟 合 模 型 的 形态 学 组 成 成 分 中 . 

图 11 展 示 了 利用 这 5 个 形态 学 成 分 的 拟 合 结 果 ， 
各 子 结构 成 分 的 模型 图 像 结果 可 见 图 12. 相 比 于 
他 解构 模式 , 添加 内 旋 臂 使 得 模型 图 像 和 观测 图 像 
最 为 一 致 ; 在 残 差 图 像 中 已 经 看 不 到 任何 明显 的 子 
结构 残留 . 这 说 明 M81 中 明显 的 子 结构 都 已 被 考虑 
进 模型 . 从 面 亮 度 径 向 轮廓 上 看 , 模型 基本 能 够 很 
好 地 预测 M81 的 面 亮度 、 椭 率 和 方位 角 的 径 向 变 
化 . 

从 表 1 中 能 够 清晰 地 看 出 添加 内 旋 臂 后 核 球 和 
盘 性 质 的 量化 变化 . 核 球 的 主要 变化 发 生 在 椭 率 和 
方位 角 上 , 分 别 变化 了 6 加 和 8%, 其 他 参数 变化 不 
明显 . 内 旋 臂 的 添加 对 盘 的 影响 非常 明显 : 添加 了 
内 旋 臂 之 后 , 盘 的 亮度 下 降 了 约 0.5 mag, 有 效 半 径 
(76) 减 小 约 38%, Sérsic 指 数 n 减 小 约 50%, 光度 占 比 
下 降 18.68%; 而 内 旋 辟 本 身 的 光度 占 比 是 17.06%. 
这 说 明了 在 模型 没有 考虑 内 旋 臂 的 情况 下 ， 
GALFIT 在 解构 时 把 一 部 分 内 旋 臂 的 光 分 配给 了 
di. 除 此 之 外 其 他 子 结构 的 参数 值 变化 不 明显 . 残 
差 图 像 的 光度 占 比 约 为 8.58%, 主要 来 自 剩 余 的 前 
景 星 和 M81 的 其 他 恒星 . 内 旋 臂 的 添加 并 不 影响 外 
旋 臂 的 拟 合 结果 . 
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图 10 ”四 成 分 ( 核 球 十 盘 十 AGN 十 外 旋 辟 ) 的 拟 合 结果 展示 . 图 像 的 格式 和 图 6 与 图 7 相同 . 


Fig.10 The fitting results of four components (bulge+ disk+ AGN+ outer spiral arms). The format of the image is the same as 
that in Fig. 6 and Fig. 7. 
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图 11 五 成 分 的 拟 合 结果 展示 . 图 像 的 格式 和 图 6 与 图 7 相同 . 


Fig.11 The fitting results of five components. The format of the image is the same as that in Fig. 6 and Fig. 7. 
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64 3$ 天 X 学 dk 2H 


1 component 
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图 12 所 有 模式 中 的 模型 和 子 结构 的 图 像 展 示 , 其 中 的 模型 图 像 是 相应 模式 中 所 有 子 结构 图 像 的 总 和 ， 


Fig.12 Images for all models and all substructures in all the decomposition types, where the model images are a sum of the 


images for all the substructures in corresponding types. 
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64 48 defe. 早 型 旋涡 星系 M81 (NGC 3031) 的 形态 学 解构 研究 2 
5 wit RBA 2E 45796.6496. 0.95%, 12.34%, 5.26% 


5.1 ”6 种 星系 解构 模式 的 主要 结果 比较 

表 1 记 录 了 所 有 的 参数 值 结果 , 为 此 我 们 可 以 
定量 对 比 不 同 星系 解构 策略 带 来 的 影响 . 对 于 模型 
的 总 星 等 , 除 单 成 分 模型 的 星 等 与 观测 图 像 的 结果 
相等 以 外 , 其 他 所 有 的 解构 策略 都 略 低估 了 0.02 星 
等 , 这 样 小 的 偏差 比 背 景 的 不 确定 性 带 来 的 误差 还 
要 低 . 也 就 是 说 所 有 的 星系 解构 策略 都 可 以 对 亮度 
进行 精准 预测 . 

我 们 的 结果 显示 实施 不 同 的 解构 模式 确实 会 
导致 星系 整体 形态 学 的 估计 产生 偏差 . 其 中 , 过 于 
简单 的 单 成 分 ( 单 S6rsic 函 数 ) 解 构 模 式 得 到 的 S6rsic 
由 数 n 为 核 球 的 值 ， 有 效 半 径 (7。) 为 盘 的 值 , 轴 比 和 
方位 角 大 约 是 核 球 和 盘 的 平均 值 . 这 说 明 如 同 M81 
这 样 的 星系 是 一 个 包含 多 个 子 结构 的 复杂 的 形态 
学 系统 , 星系 的 整体 性 质 是 不 同 子 结构 各 自 的 性 质 
的 综合 . 

除了 单 成 分 解构 模式 , 其 他 5 种 形态 学 解构 对 
核 球 的 形态 学 参数 值 都 没有 十 分 明显 的 变化 . 以 每 
个 参数 的 最 小 值 为 基准 , 核 球 的 光度 占 比 、 星 等 、 
了 效 半径 、Sérsic 指 数 、 轴 比 和 方位 角 的 最 大 变化 
EE 31E6.9096. 0.7796. 9.4096. 6.3696. 7.14% 
和 13.10%, 其 中 除了 方位 角 以 外 , 核 球 的 参数 变化 
比 背 景 的 不 确定 性 带 来 的 误差 更 低 . 对 于 方位 角 ， 
核 球 上 + 盘 的 双 成 分 解构 和 核 球 + 盘 +AGN 的 三 成 分 
解构 得 到 的 结果 一 致 ; 而 核 球 + 盘 +AGN 十 外 旋 臂 
的 四 成 分 解构 和 核 球 + 盘 +AGN+ 外 旋 臂 二 内 旋 辟 
的 五 成 分 解构 分 别 升 高 8.41% 和 降低 4.41%. 这 个 结 


N 


= 


2.08% LA 10.59%. 而 在 添加 了 内 旋 臂 的 五 成 分 
解构 结果 中 , 盘 的 上 述 参 数 分 别 变化 了 23.15%、 
7.30%、24.89%、85.71%、4.25% 以 及 26.06%. 因此 
可 以 得 出 , 除了 内 旋 臂 , 其 他 子 结构 是 否 被 单独 拟 
合并 不 能 明显 影响 盘 的 参数 值 . 对 于 方位 角 , 在 不 
同 解构 模式 中 , 盘 的 变化 比 核 球 大 , 这 是 因为 盘 与 
旋 辟 直接 重合 . 内 旋 臂 是 否 单独 拟 合 对 盘 的 各 参数 
影响 最 大 , 这 是 因为 如 果 模 型 中 不 单独 考虑 内 旋 
E, 那么 内 旋 臂 的 光 会 被 分 配 于 盘 . 
上 述 的 所 有 模式 的 实验 都 表明 AGN 对 星系 整 
体 以 及 各 子 结构 的 拟 合 影响 不 大 . 在 这 6 种 星系 解 
构 模 式 当 中 , AGN 的 星 等 变化 (以 最 小 值 为 基准 ) 都 
在 14.33% 以 内 , 光度 占 比 的 变化 也 在 0.52% 以 内 . 这 
说 明 对 AGN 亮 度 的 估计 也 同样 不 易 受 到 其 他 子 结 
构 的 影响 . 但 是 AGN 在 1 个 Skrsic 函 数 +TAGN 的 双 成 
分 解构 模式 结果 中 比 其 他 模式 结果 低 两 个 星 等 , 这 
是 因为 没有 添加 核 球 成 分 而 导致 对 AGN 亮 度 的 高 
fh. 


5.2 不同 星系 解构 模式 的 适用 性 

通过 上 述 多 种 解构 模式 的 实验 , 我 们 发 现 对 于 
复杂 度 越 高 的 解构 模式 GALFIT 运 行 的 时 间 越 长 , 特 
别 是 有 旋 臂 建 模 的 坐标 旋转 函数 以 及 对 截断 建 模 
的 截断 函数 参与 的 解构 , 由 于 参数 空间 大 , 因此 每 
次 拟 合 都 需要 花费 比 简 单 模型 更 久 的 时 间 , 对 计算 
机 配置 的 要 求 也 要 更 高 . 对 旋 臂 或 截断 的 建 模 需 要 
人 为 调节 参数 , 这 种 模式 对 大 样本 星系 的 研究 显然 
是 有 局 限 性 的 . 此 外 , 我 们 的 结果 显示 , 添加 了 外 旋 


FE ME 


果 表 明 对 核 球 方位 角 的 估计 易 受 旋 臂 的 影响 ; 即便 


臂 和 内 旋 臂 之 后 , 核 球 和 盘 的 参数 也 会 发 生 不 同 程 


外 旋 臂 不 与 核 球 重合, 也 会 导致 核 球 方位 角 的 估计 
偏差 大 于 背景 不 确定 性 . 

相 比 于 核 球 , 盘 结构 在 不 同 的 解构 模式 中 性 
质变 化 比较 大 . 以 各 参数 的 最 小 值 为 基准 , 盘 的 
光度 占 比 、 星 等 、 有 效 半 径 、Sérsic 指 数 、 轴 比 、 
方位 角 的 最 大 变化 幅度 分 别 为 64.52%、7.30%、 
38.3096. 93.5496. 4.269604 7€ 26.0696. 核 球 + 盘 的 
双 成 分 和 核 球 + 盘 +AGN 的 三 成 分 解构 结果 基本 
一 致 . 相 比 于 核 球 + 盘 的 双 成 分 或 核 球 + 盘 +AGN 
的 三 成 分 解构 , 添加 了 外 旋 臂 的 四 成 分 解构 使 上 
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度 的 变化 . 因此 本 工作 可 以 为 其 他 星系 形态 学 解构 
研究 提供 一 些 经 验 参 考 . 

在 核 球 + 盘 的 双 成 分 解构 和 核 球 + 盘 +AGN 的 
三 成 分 解构 中 , 核 球 和 盘 的 结果 相差 不 大 ; 从 我 们 
的 经 验 来 看 , 核 球 + 盘 的 双 成 分 解构 和 核 球 二 盘 十 
AGN 的 三 成 分 解构 结果 均 非常 稳定 , 不 易 受 初始 
参数 的 影响 ; 即使 初始 参数 的 估计 并 不 合理 , 这 两 
种 解构 模式 也 会 得 到 一 致 合 理 的 结果 ; 这 两 种 解构 
的 运行 过 程 耗 时 最 短 ( 除 了 过 于 简单 的 1 个 SErsic 函 
数 的 单 成 分 解构 和 1 个 Srsic 函 数 +TAGN 的 双 成 分 
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解构 ). 通过 比较 这 两 种 解构 模式 的 结果 可 以 发 现 ， 
是 否 添 加 PSF 函 数 对 核 球 盘 的 拟 合影 响 不 大 . 不 
过 我 们 建议 添加 PSF 函 数 来 拟 合 星系 核心 的 结构 . 
因为 虽然 对 于 M81 来 说 AGN 亮 度 不 高 , 但 是 对 于 那 
些 AGN 比 较 亮 (或 PSF 存 在 比较 明显 外 晕 ) 的 星系 
图 像 来 说 , PSF 的 添加 将 提高 星系 解构 的 准确 性 ， 
而 添加 PSF 函 数 对 整个 解构 的 耗 时 影响 非常 小 . 因 
此 , 我 们 建议 核 球 + 盘 +AGN 的 三 成 分 拟 合 是 最 适 
于 大 样本 星系 形态 学 解构 的 模式 . 

在 我 们 的 工作 中 , 核 球 + 盘 +AGN 
四 成 分 解构 与 核 球 + 盘 +AGN 二 外 旋 辟 
五 成 分 解构 对 于 还 原 星系 整体 的 观测 特征 来 说 最 
为 准确 , 不 过 这 两 种 解构 过 程 耗 时 最 长 , 因此 比较 
| 个 源 星系 进行 深入 细致 研究 的 工作 , 比如 


外 旋 臂 的 
Vi Jie ^ Bh] 


适用 于 对 
单独 测量 星系 旋 臂 的 参数 等 . 有 研究 认为 位 于 核 球 
以 外 的 子 结构 (如 旋 辟 等 ) 对 核 球 的 形态 学 参数 的 拟 
合 结果 影响 不 大 ( 相 比 那些 位 于 核 球 附近 子 结构 如 
棒 等 )523. 我 们 的 结果 部 分 支持 这 种 结论 , 但 也 有 不 


同 : 相同 点 是 我 们 的 核 球 + 盘 +AGN+ 外 旋 臂 的 四 
成 分 解构 和 核 球 + 盘 +AGN 二 外 旋 臂 上 + 内 旋 臂 的 五 
成 分 解构 结果 显示 外 旋 臂 的 独立 拟 合 与 否 确实 对 
核 球 和 盘 的 影响 都 很 小 , 而 在 五 成 分 解构 结果 中 核 
球 附 近 内 旋 臂 独立 拟 合 与 否 对 核 球 的 影响 也 很 小 
但 对 盘 的 影响 却 很 大 . 因此 我 们 建议 , 对 于 那些 位 
于 核 球 附近 的 子 结构 需要 酌情 考虑 , 是 否 将 它们 独 
立 拟 合 可 能 跟 核 球 以 外 的 子 结构 一 样 也 不 影响 对 
核 球 的 拟 合 , 但 是 会 影响 盘 结 构 的 拟 合 . 


5.3 ”AGN 的 存在 对 星系 解构 的 影响 

AGN 作 为 星系 中 心 的 子 结构 , 其 亮度 可 能 会 
对 寄主 星系 的 测量 造成 一 定 的 影响 . 本 工作 很 重要 
的 一 部 分 内 容 就 是 探究 AGN 在 不 同 星系 解构 模式 
中 的 影响 . 为 此 , 我 们 实际 还 额外 做 了 核 球 上 + 稳 十 外 
旋 臂 + 内 旋 臂 四 成 分 解构 实验 . 图 13 记 录 了 这 个 实 
验 的 图 像 结果 , 其 中 各 个 成 分 的 模型 图 像 见 图 12. 
表 1 中 也 记录 这 个 实验 的 参数 值 . 


AJ 


(1) Observation 


(4) Surface brightness profile 


(3) Residual 
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13 ”无 AGN 的 四 成 分 ( 核 球 十 盘 十 外 旋 辟 十 内 旋 辟 ) 的 拟 合 结果 居 


图 


RINK SAI 


AJ. 7 相同 . 


Fig.13 The fitting results of four components (bulge 十 disk + outer spiral arms + inner spiral arms). The format of the image 


is the same as that in Fig. 6 and Fig. 7. 
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型 旋涡 星系 M81 (NGC 3031) 的 形态 学 解构 


九 


与 核 球 + 盘 +AGN+ 外 旋 臂 + 内 旋 臂 的 五 成 分 
解构 相 比 我 们 发 现 , AGN 成 分 的 独立 拟 合 与 否 对 
各 子 结构 成 分 的 参数 影响 都 非常 小 . 最 大 变化 是 核 
球 的 Srsic 指 数 , 添加 了 AGN 之 后 核 球 的 S6érsic 指 数 
下 降 约 4.70%; 其 他 的 参数 在 添加 AGN 之 后 的 变化 
都 在 1% 以 内 . 所 有 参数 的 变化 都 比 背 景 不 确定 性 
带 来 的 误差 小 很 多 . 模型 图 像 和 残 差 图 像 的 对 比 看 
不 出 明显 差别 , 径 向 面 亮 度 轮廓 同样 没有 明显 变 
化 . 残 差 图 像 光 度 占 比 约 8.61%, 和 五 成 分 的 残 差 图 


5.5 星系 解构 过 程 中 的 常见 问题 

在 最 后 的 讨论 部 分 , 我 们 为 读者 提供 一 些 在 
本 工作 中 汲取 的 经 验 . 首先 , 我 们 的 结果 建议 在 开 
始 拟 合 之 前 要 了 解 GALFIT 软 件 中 图 像 和 参数 的 设 
置 , 例如 通常 望远镜 数据 官网 提供 的 权重 (Weight) 
图 像 、 误 差 图 像 或 噪声 图 像 可 能 和 GALFIT 的 要 求 
不 同 , 图 像 输 入 的 格式 、 单 位 与 要 求 不 一 致 会 得 到 
错误 的 结果 ; 其 次 , 我 们 建议 拟 合 窗口 (或 视 场 ) 设 
置 得 大 一 些 , 因为 GALFIT 构 建 出 的 模型 的 质量 对 背 


像 类 似 , 这 些 光 主要 来 自 剩 余 的 前 景 星 和 星系 中 的 
其 他 恒星 . 

我 们 的 实验 覆盖 了 从 简单 到 复杂 的 多 种 模 
来 讨论 AGN 对 不 同 解构 模式 的 影响 , 得 到 的 结果 


b 


RE 


£H ZR 
是 添加 AGN 与 否 对 各 子 结构 参数 影响 均 保持 在 
约 1% 以 内 , 比 背 景 不 确定 性 带 来 的 误差 小 . 尽管 如 
此 , 我 们 仍然 建议 独立 拟 合 AGN (或 星系 核心 ), Al 
为 M81 的 AGN 是 低 光度 的 , 所 以 它 对 解构 几乎 不 产 
生 明 显 的 影响 . 然而 , 对 于 高 光度 AGN, 预期 影响 
还 是 很 可 观 ; 即便 未 知 AGN 的 光度 或 是 否 存 在 , 我 
们 也 推荐 独立 拟 合 星系 核心 . 


5.4 ”更 加 复杂 的 解构 模式 

在 本 工作 中 , 我 们 采用 了 不 同 复杂 程度 的 6 种 
系 解构 模式 , 但 是 我 们 并 不 认为 模型 包含 的 形态 
学 成 分 越 多 越 好 . 本 工作 最 复杂 的 模型 是 五 成 分 子 
结构 模型 . 在 五 成 分 解构 结果 的 残 差 图 像 中 不 能 明 
显 看 见 还 有 其 他 形态 学 成 分 , 但 并 不 能 由 此 下 结 
论说 M81 确 实 由 这 5 个 子 结构 成 分 组 成 . 有 研究 表 
明 M81 的 内 、 外 旋 臂 本 是 一 体 , 同属 于 一 种 旋涡 状 
子 结构 69-60. 不 过 在 M81 的 紫外 图 像 上 看 不 见 内 
旋 臂 ,这 至 少 说 明 内 、 外 旋 臂 有 不 同 的 星 族 成 分 . 
因此 , 在 我 们 的 工作 中 虽然 区 分 了 内 、 外 旋 辟 , 但 
不 能 证 明 它们 是 独立 的 子 结构 . 

确实 , 添加 更 多 的 成 分 可 能 会 得 到 更 加 理想 的 

拟 合 结果 , 但 是 这 些 子 结构 的 动力 学 归属 并 非 是 形 
态 学 解构 可 以 回答 的 问题 . 此 外 , 由 于 一 些 子 结构 
的 亮度 有 限 , 对 星系 的 主要 成 分 (如 核 球 和 盘 ) 几 乎 
没有 影响 , 因此 我 们 建议 , 只 要 研究 目标 不 是 为 了 
专门 研究 这 些 暗 弱 的 子 结构 , 在 星系 解构 过 程 中 不 
必 独 立 拟 合 这 些 低 亮度 的 潜在 子 结构 . 
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景 敏感 (不 过 需要 注意 大 窗口 会 降低 拟 合 的 速度 ， 
因此 需要 权衡 拟 合 窗口 的 大 小 ); 再 次 , 需要 常用 
GALFIT 的 模型 输出 功能 , 特别 是 复杂 解构 的 时 候 ， 
提前 人 为 构建 一 个 与 目标 相像 的 模型 能 极 大 地 提 
高 解构 效率 . 最 后 , 可 能 是 GALFIT 的 一 个 潜 规 则 , 输 
出 的 文件 名 必须 以 字母 开头 而 不 是 数字 , 否则 拟 合 
无 法 运行 . 希望 这 些 经 验 能 对 读者 有 所 帮助 . 


总 结 


IA ZH 


本 工作 是 一 项 基于 Spitzer-IRAC 4.5 hm 图 像 ， 
利用 星系 解构 软件 GALFIT 对 近邻 早 型 旋涡 星系 
M81 进 行 形态 学 解构 的 研究 , 目的 是 量化 M81 中 
子 结构 的 形态 参数 并 探究 不 同 解构 模式 对 核 球 和 
盘 等 主要 子 结构 拟 合 的 影响 . 我 们 采用 了 由 简单 到 
复杂 的 6 种 形态 学 解构 模式 , 结果 显示 , 星系 整体 的 
形态 ( 即 不 独立 拟 合 任何 子 结构 ) 是 星系 中 存在 的 所 
有 子 结 构 形 态 的 综合 , 只 有 进行 形态 学 解构 , 才能 
全 面 、 完 善 地 认识 星系 的 形态 以 及 其 他 相关 性 质 . 
本 工作 最 主要 的 结论 如 下 : 
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e M81 的 核 球 是 一 个 Sérsic 指 数 约 为 ”有效 
半径 约 1.3 kpce、 轴 比 约 为 0.7、 方 位 角 约 为 
一 34*、 视 星 等 约 为 7.7 mag 的 经 典 核 球 , 其 光 
度 占 比 约 为 40%. 核 球 的 这 些 参 数值 在 所 有 的 
解构 模式 中 均 保 持 稳 定 , 不 易 被 其 他 子 结构 


影响 . 


对 于 盘 结 构 来 说 , 不 同形 态 学 解构 模式 得 到 
了 不 同 的 结果 . 除了 五 成 分 拟 合 , 其 他 解构 模 
式 均 得 到 M81 的 稻 的 Sérsic 指 数 约 为 1.2、 有 
效 半径 约 为 3.8 kpc、 轴 比 约 为 0.47、 方 位 角 
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约 为 -26*、 视 星 等 约 为 7.4 mag, 其 光度 占 比 
71735696, 它 是 M81 在 4.5 hm 图 像 中 最 亮 的 子 
结构 . 但 是 将 内 旋 臂 独立 拟 合 后 , 盘 结构 变 
得 更 加 平坦 , 这 时 Skrsic 指 数 约 为 0.6、 有 效 
半径 约 为 5 kpe, 光度 占 比 约 为 35%, 是 M81 在 


4.5 pm 波段 第 2 亮 的 子 结构 . 


旋 臂 的 形态 学 参数 相对 稳定 . 外 旋 臂 和 内 旋 
臂 的 星 等 分 别 约 为 9.8 mag 和 8.7 mag, 光度 占 
EE 4 310 25 7369611796, 这 两 种 子 结构 都 不 易 
受到 其 他 子 结构 的 影响 . 


IRA 

AGN 在 所 有 形态 学 解构 模式 中 的 结果 都 相 
对 稳定 . 除了 单 S6rsic 函 数 二 PSF 的 双 成 分 拟 
合 , AGN 的 星 等 均 约 为 14 mag, 光度 占 比 约 为 
0.196. 而 在 单 SGrsic 函 数 +PSE 的 双 成 分 拟 合 
结果 中 AGN 的 亮度 更 大 (其 星 等 约 12 mag, 光 
度 占 比 约 为 0.64%), 是 由 于 没有 核 球 的 独立 
拟 合 而 导致 亮度 的 高 估 ( 包 含 部 分 核 球 光度 ). 


通过 对 不 同形 态 学 解构 模式 的 结果 进行 比较 ， 
我 们 建议 解构 模式 的 选择 取决 于 具体 研究 的 目标 . 
大 样本 研究 适合 采用 核 球 + 盘 +AGN 的 三 成 分 拟 
合 , 而 对 于 个 源 星 系 的 研究 则 适合 更 多 成 分 (比如 
添加 旋 臂 等 子 结构 ) 的 独立 拟 合 . 此 外 , 我 们 发 现 
旋 臂 (特别 是 内 旋 臂 , 如 有 的 话 ) 对 盘 的 影响 非常 大 ， 
而 对 核 球 的 影响 很 小 , 因此 我 们 建议 如 果 想 要 精确 
测量 盘 的 形态 参数 则 需要 考虑 旋 臂 的 影响 . 

这 项 工作 是 我 们 对 M81 一 系列 研究 的 开始 . 接 
下 来 我 们 会 进行 M81 子 结构 的 多 波段 形态 学 研究 ， 
量化 其 形态 学 K- 改 正 (Morphological K-Correc- 
tion), 并且 运 用 SED 拟 合 来 探究 M81 星 系 整体 及 
其 子 结构 的 星 族 成 分 的 空间 分 布 (二 维 性 质 ). 在 此 
基础 上 , 我 们 将 尝试 利用 数值 模拟 演绎 M81 及 其 子 
结构 的 形成 和 演化 . 


致谢 我 们 感谢 审 稿 人 对 这 篇 文章 细致 的 审阅 并 提 
出 具有 建设 性 意义 的 意见 , 这 些 意见 让 本 文 收获 非 
常 大 的 改进 . 我 们 感谢 广州 大 学 天 体 物理 中 心 全 体 
老师 对 本 工作 的 建议 .在 工作 中 我 们 利用 了 NASA/ 
IPAC 河 外 星系 数据 库 (NED) 以 及 SAO/NASA 天 体 
物理 数据 系统 (ADS). 
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Asstract This paper presents a morphological investigation of the early-type spiral galaxy M81 (NGC 
3031) through decomposition by fitting radial profiles of surface brightness using the software GALFIT, 
aimed at exploring structural components of M81 and quantifying their morphology. In this work, we 
adopt 6 types of decompositions with different numbers of morphological components, among which the 
most complicated one contains 5 components such as a bulge, a disk, a pair of outer spiral arms, a pair 
of inner spiral arms, and a galactic nucleus. The results show that, M81 hosts a classical bulge with 
the Sérsic index ~ 5.0; the morphology and the luminosity for the bulge are almost constant in the 
different decomposition types. The disk of M81 has the Sérsic index ~ 1.2 but the morphology and the 
luminosity are found sensitive to the inclusion of the inner spiral arms in the decomposition or not. The 
results of this work indicate that the combination of individual substructures has a considerable impact 
on the morphology of the galaxy as a mixture. On the basis of the results, the usability of the different 
types of decomposition is suggested in this work. The three-component decomposition, i.e., bulge + disk 
+ nucleus, is applicable to statistics of large samples of galaxies; more complicated decomposition with 
spiral arms taken into account is suitable for precise measurements of individual galaxies in small samples. 
This work is based on the single-band Spitzer-The Infrared Array Camera (IRAC) 4.5 um image. In the 
future, we will carry out multi-wavelength decomposition, for the purpose of investigating spectral energy 
distributions and stellar population properties for the galactic substructures, and thereby deduce their 
formation history and evolution processes. 


Key words galaxies: individual: M81, bulges, galaxies: fundamental parameters: morphology, galaxies: 
structure, galaxies: spiral 
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